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First measurement of the helicity asymmetry E in 77 photoproduction on the proton 
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Abstract 

Results are presented for the first measurement of the double-polarization helicity asymmetry E for the rj photoproduction 
reaction yp —► r/p. Data were obtained using the FROzen Spin Target (FROST) with the CLAS spectrometer in Hall B 
at Jefferson Lab, covering a range of center-of-mass energy W from threshold to 2.15 GeV and a large range in center- 
of-mass polar angle. As an initial application of these data, the results have been incorporated into the Julich-Bonn 
model to examine the case for the existence of a narrow N* resonance between 1.66 and 1.70 GeV. The addition of 
these data to the world database results in marked changes in the predictions for the E observable from that model. 
Further comparison with several theoretical approaches indicates these data will significantly enhance our understanding 
of nucleon resonances. 

Keywords: eta photoproduction; polarization observable; helicity assymmetry 


1. Introduction 

Much activity is being devoted to establishing the de¬ 
tails of the excitation spectrum of the nucleon in order to 
deepen our understanding of that fundamental strongly- 
interacting three-quark system. Due to the broad widths 
of the nucleon excitations (of the order of 100-300 MeV), 
the states overlap in the mass spectrum. Thus, disentan¬ 
gling the individual states to identify their exact masses 
and quantum numbers has been quite difficult. While 
some resonances are well established, fewer states have 
been observed than most constituent quark models and 
Lattice QCD calculations predict |1]. An additional com¬ 
plexity arises because, beyond resonance states with typi¬ 
cal widths, approaches based on chiral quark solitons also 
predict states with far narrower widths than do constituent 
quark models, including, for example, an N^ + state with 
a width of 40 MeV or less [H OH 010 [5] at about 1.7 GeV; 
this particular state may have been observed in r] photo¬ 
production on the neutron mum- 

Since differential cross section measurements alone are 
insufficient to locate the underlying resonance states or 
determine their properties, attention has turned to polar¬ 
ization observables. Polarization observables involve inter¬ 
ferences between sets of amplitudes, so their measurement 
can provide stringent tests for predictions of the photo¬ 
production process and help sort out ambiguities in the 
theoretical description of the reaction in terms the reso¬ 
nances involved. One such polarization observable is the 
helicity asymmetry E in pseudoscalar meson photopro¬ 
duction, which is the normalized difference in photopro¬ 
duction yield when spins of the incident photon and a 
longitudinally-polarized target are parallel and anti-parallel. 


Formally, this observable is defined as a modulation of the 
center-of-mass differential cross section da/d^lo through 
the relation 


da 

dfl 


da 

dflo 


(1 - PjP^E) , 


( 1 ) 


where Pj specifies the degree of longitudinal target polar¬ 
ization and Po is the circular polarization fraction of the 
incident photon beam. This asymmetry is generally ex¬ 
pressed as a function of the center-of-mass energy W and 
the polar angle of the produced meson in the center-of- 
mass frame cos 9 cm . 

Pion photoproduction studies have contributed greatly 
to the knowledge of the nucleon resonance spectrum. Re¬ 
cent measurements of spin observables in pion photopro¬ 
duction nnunjEunam], have illustrated the power of 
polarization observables to clarify that spectrum. Even 
so, many ambiguities still exist and many predicted states 
remain unobserved. Though pion photoproduction offers 
a larger cross section, the photoproduction of p mesons ex¬ 
hibits the interesting feature that the process excludes con¬ 
tributions from resonances with isospin I = 3/2, thereby 
isolating the N* (7=1/2) states. The t) photoproduction 
process on the proton thus acts as an “isospin filter” for the 
nucleon resonance spectrum, resulting in a useful tool for 
disentangling the different states, and is especially impor¬ 
tant in finding and investigating states that do not couple 
strongly to pions. 


2. The Experiment 

The measurements reported here are an integral part of 
a program at Jefferson Lab to achieve a “complete” exper- 
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iment for the 77 pliotoproduction process, whereby all the 
helicity amplitudes are determined for photoproduction of 
that pseudoscalar meson. The program began with mea¬ 
surements of the unpolarized differential cross section 
mm using the large solid angle CEBAF Large Accep¬ 
tance Spectrometer (CLAS) [17) . and the bremsstrahlung 
photon tagger housed in Jefferson Lab Hall B [15] , For 
the measurements reported here, circularly polarized pho¬ 
ton beams were produced by polarization transfer from the 
polarized electron beam of the CEBAF accelerator, which 
was incident on an amorphous radiator of the photon tag¬ 
ger. 

The target nucleons for the pliotoproduction process 
were free protons in frozen butanol (C 4 H 9 OH) beads in¬ 
side a 50-mm-long target cup. [19]. The protons of the 
hydrogen atoms in this material were dynamically polar¬ 
ized along the photon beam direction. The longitudinal 
target polarization PT was determined with nuclear mag¬ 
netic resonance measurements, and averaged 82 ± 5%. To 
minimize systematic uncertainties, the orientation of the 
target polarization direction was flipped every few days of 
data-taking between being aligned and anti-aligned with 
the direction of the incoming photon beam. The helicity 
of the beam was flipped at a rate of 30 Hz. 

Final state particles resulting from pliotoproduction 
were detected using CLAS, a set of six identical charged 
particle detectors installed in a toroidal magnetic field. 
The principal CLAS subsystems required for this study 
were the drift chamber systems for tracking charged parti¬ 
cles J 2 D 1 , a scintillator-based time-of-flight system ED, and 
a start counter array which determined when charged par¬ 
ticles passed from the target into the detection region [ 22 ] . 
The energy and polarization information for incident pho¬ 
tons was provided by the photon tagger. 

3. Analysis 

To determine the helicity asymmetry £ in a discrete 
event counting experiment, Eq. ([T]) is inverted to form the 
asymmetry 


1 /A+-iV_ 

~~\PT\\P°\ \N+ + N- 

where the detector acceptance cancels. The cross sections 
are replaced by N + and _ZV_, which are the number of 77 
mesons counted in beam-target helicity aligned and anti¬ 
aligned settings, respectively. The background from 11011-77 
final states and those from events arising from the unpo¬ 
larized nucleons within the butanol are subtracted before 
forming this asymmetry. 

Determination of the E observable requires knowledge 
of the degree of polarization for both the photon beam 
and the target proton. The photon beam polarization is 
calculated from the incident photon energy E 1 relative to 
the bremsstrahlung endpoint (E = E 1 /E e -) via the ex¬ 



pression 


PI = P, 


4E- E 2 
4 — 4E + 3E 2 ’ 


(3) 


where P e is the polarization of the electron beam incident 
on the amorphous radiator within the photon tagger [23]; 
P e was measured with the Hall B Mpller polarimeter dur¬ 
ing the experiment to be 0.84 ± 0.01. 

Events in the detector were reconstructed in the follow¬ 
ing manner. Individual charged tracks were reconstructed 
in the CLAS drift chambers and matched to hits in the 
time-of-flight (TOF) and start counter paddles. The parti¬ 
cle identity was determined by combining the information 
on the momentum of the particle, which was determined 
by the drift chambers from the curvature of the particle 
trajectory in the magnetic field, and on the speed of the 
particle (f3) as determined from the timing information 
provided by the tagger, start counter, and TOF systems. 
Charged tracks that could not be reconstructed by all of 
these detectors were rejected. A track was assumed to 
have the particle identity that allowed the closest match 
between the 4-momentum-computed /3 and the measured 
value of /?. An additional requirement that the measured 
/3 was within ±0.04 of the expected value was imposed 
on pion candidates, significantly suppressing the electron 
background. Once the particle identity was established, 
a correction due to energy loss in the target and detector 
materials was performed, with the 4-vector values adjusted 
accordingly. The tracks and the event as a whole were as¬ 
sociated to beam photons based on consistency with the 
projected vertex timing. To avoid ambiguity, only events 
with particles matching exactly one beam photon were 
kept. 

The CLAS detector is primarily a charged particle spec¬ 
trometer, with electromagnetic calorimetry confined to a 
narrow angular range. Thus, ~ 94% of the signal in this 
analysis relied on missing mass reconstruction of the neu¬ 
tral 77 from the measured kinematical information of the 
proton recoiling into the CLAS (the detection of which 
was required), with the remainder of events having one or 
both charged pions from the decay 77 —> detected. 

Events with a single detected charged pion were required 
to have a missing mass squared greater than 0.06 GeV 2 /c 4 , 
which is the onset of the remaining two-pion phase space. 
Events with both a 7 r + and n~ detected were required to 
have the remaining missing mass squared close to that of 
the 7T° within the detector resolution: 0.008—0.028 GeV 2 /c 4 

The 77 photoproduction data were analyzed to extract 
the helicity asymmetry E in 50MeV-wide center-of-mass 
energy W bins and 0.2-wide center-of-mass production 77 
polar angle (cos0 cm ) bins. Binning in W begins near the 77 
threshold at 1.5 GeV. These bin widths were chosen to bal¬ 
ance between minimizing statistical uncertainties for the 
extraction while achieving the best energy resolution for 
the resonance spectrum and most thorough knowledge of 
the polar distribution of the resonance decay. The anal¬ 
ysis procedures described below were performed for each 
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Figure 1: Analysis example for the kinematic bin (1650 < W < 1700MeV, —0.2 < cos0 cm < 0.0). Left panel: Background fit to the missing 
mass spectra for the two helicities (higher amplitude N— is green). Middle panel: Background subtraction and net ri yield. Right panel: 
Yield difference, with fit to sidebands to determine overall asymmetry offset. (Color online) 


kinematic bin separately. 

To distinguish a photoproduced 77 from the background, 
fits were performed to the invariant mass spectra with 
models of the signal and background included, as shown 
in Fig. |T] The integral of the fit shape of the background 
and the uncertainty of this integral from the error matrix 
estimated the background contribution. Since the detec¬ 
tor resolution dominated the shape of the 77 enhancement, 
the signal was modeled as a Gaussian. Polynomials were 
used to model the background, with the order of the poly¬ 
nomial increased progressively with every fit iteration up 
to fifth order as long as the fit improvement was statis¬ 
tically significant. Specifically, improvement in the confi¬ 
dence level beyond 0.5 was considered not significant. For 
the two W bins near the 77 threshold (W < 1.6 GeV), the 
step function-like drop-off in photoproduced system phase 
space required a different approach. For those bins, the 
error function erf was used in addition to the polynomial, 
with its amplitude and transition width as free parameters. 

A single fit was performed to the spectra of both beam 
helicities, with a common model of the background shape 
and common position and width of the 77 enhancement; an 
example is shown in the left panel of Fig. [T] The unpo¬ 
larized background essentially cancels out in the difference 
of the yield in the missing mass spectra of the two helic¬ 
ities ( N + and 1V_) seen in the numerator of Eq. ([ 2 ]); an 
example is given in the middle panel of Fig. [l] The small 
remaining overall vertical offset seen for some kinematic 
bins may be due to asymmetries in the broad polarized 
background, such as the asymmetry which might exist for 
the 7 r + 7 r - final state. These offsets were determined with 
a fit to the sidebands (as seen in the right panel of Fig. [l]), 
defined as ±3er from the peak center, where the center 
and a values were derived from the previously performed 
fit of the 77 enhancement. The normalized asymmetry was 
thus calculated from this corrected difference of helicities 
divided by the overall 77 yield determined with the back¬ 


ground subtraction described above. 

A separate study of photoproduction on a pure carbon 
target showed no evidence of peaking in the 77 mass region. 
Thus, no correction for the heavier nuclei in the target was 
required beyond the smooth background fitting described. 
Helicity asymmetry extraction was not performed when 
the background exhibited an extremum under the 77 peak, 
(z.e. within ±lcr of the peak centroid) to avoid serious 
ambiguities between the signal and the background shape. 
Additionally, analysis in a kinematic bin was abandoned 
when the total 77 yield uncertainty was greater than 30%. 

Missing mass energy resolution for the 77 with CLAS is 
a smooth function of the kinematic space explored here. 
Therefore, the peak widths seen in the initial independent 
analyses of the individual kinematic bins were compared 
and a smooth function of the peak width across the kine¬ 
matic space was extracted. Yield extractions were then 
repeated using these constraints on the peak width to en¬ 
force consistency with the detector resolution. 

Statistical uncertainties dominated the systematic un¬ 
certainties in all analyzed bins, and are shown combined 
in the presented results. The systematic uncertainties in¬ 
clude the target polarization Pj uncertainty 6 . 1 %) and 
photon beam polarization Po uncertainty (3.1%). 

4. Results and discussion 

The results for the helicity asymmetry E are shown in 
Figs. [2] and [3] for 1.5 < W < 2.1 GeV. At threshold, the 
E observable is close to unity due to the dominance of the 
1V(1535)| resonance [2i] . and the results reported here 
are consistent within uncertainties with this expectation. 
As W increases, the presence of other resonances and the 
interferences of the various amplitudes related to those 
resonances generate a TV-dependent structure in E, which 
models of the production process attempt to describe. As 
examples of such models, shown in Fig. [2] are predictions 
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from phenomenological fits by SAID [55], Jiilich-Bonn [5B] 
and ANL-Osaka m- 

The figure also shows a new fit with the Jiilich-Bonn 
dynamical coupled-channel approach incorporating the data 
reported here. In that framework, the hadronic scatter¬ 
ing amplitude is constructed with a potential, generated 
from an effective SU(3) Lagrangian, using time-ordered 
perturbation theory, and the amplitude is iterated in a 
Lippmann-Schwinger equation such that unitarity and an- 
alyticity are automatically respected. 

This new fit also simultaneously incorporated the world 
databases for the pion-induced production of r)N, K A and 
KYi final states J2S] and the partial-wave solution of the 
GW/DAC group [25] for elastic nN scattering. It also in¬ 
cludes the world data bases of pion and rj photoproduction 
off the proton up to IT ~ 2.3 GeV [30l 126] , in particular 
the recent MAMI results on T and F in ry photoproduc¬ 
tion [51] . 

This new fit also simultaneously incorporated measure¬ 
ments for pion-induced reactions and for pion- and ry pho¬ 
toproduction off the proton. The published Jiilich-Bonn 
predictions for E, as well as this new fit, include the recent 
MAMI results on T and F in r/ photoproduction 

In order to achieve a good fit result, all parameters tied 
to the resonance states and to the photon interaction had 
to be modified from the values reported in Ref. !26j. The 
inclusion of these new E data also resulted in significant 
changes in the extracted resonance pole positions. For ex¬ 
ample, with these new E data, the _/V(1710)l/2 + resonance 
becomes 45 MeV heavier and 20 MeV wider compared to 
Ref. [251, with a 40% smaller branching ratio into the r/N 
channel. Helicity couplings for the high-spin iV(2190)7/2” 
and IV(2250)9/2“ resonances, whose properties are diffi¬ 
cult to determine in general, change their pole positions by 
up to 80 MeV in the real and 100 MeV in their imaginary 
parts due to these new E data. 

That the E data reported here have such a large im¬ 
pact on the resonance parameters might appear surprising. 
Since the different spin observables have differing combi¬ 
nations of amplitudes, the various observables will have 
differing degrees of impact on reducing the uncertainties 
of the parameter values extracted from a fit to the exper¬ 
imental data. In the present case, those parameters are 
the fundamental electromagnetic properties of resonances, 
their helicity couplings at the pole. 

To study the variations in the statistical impacts on the 
parameters of the new Jiilich fit that arise through the use 
of measurements of the different spin observables, we have 
studied the condition number n for the covariance matrix 
found in fitting the observed data. The condition num¬ 
ber is a standard test for diagnosing multicollinearity 
and, hence, the non-orthogonality of the model parame¬ 
ters on which the fit is based [32]. The condition number 
k is defined as A max /A m j n , that is, the ratio of the largest 
and smallest eigenvalues. Geometrically, yjk determines 
the ratio of the longest half-axis of the statistical uncer¬ 
tainty ellipsoid divided by the shortest one. A large k (say, 


greater than 100 ) is a sign of moderate to strong mul- 

ticollinearity, i.e. a very elongated statistical uncertainty 
ellipsoid; larger values of k thus connote greater uncer¬ 
tainty in the corresponding helicity couplings determined 
from the fit. 

For the present E data, we found n ~ 50, while for 
the other spin observables (and even the differential cross 
section), k ranged from 50 to 400. Thus, in terms of min¬ 
imizing the uncertainties in the extracted parameters, the 
E observable measurements reported here indeed turn out 
to be particularly impactful. This underscores that the 
observable E in ry photoproduction is especially suited to 
disentangle electromagnetic resonance properties. With 
relatively few data points, this measurement offers a larger 
impact on the baryon spectrum, helicity couplings, and 
even hadronic decay parameters than might be expected. 

Turning next to the putative N^ + resonance near W ~ 
1.7 GeV, Fig. [3]shows our results for the observable E us¬ 
ing finer W bins (of 20MeV width). Coarser, 0.4-wide 
binning in cos 9 cm was used to compensate for the nar¬ 
row energy binning. A fit to this re-binned data using the 
Julich-Bonn formalism found that the structure observed 
at ~ 1.7 GeV for the cos0 cm bin centered at 0.2 is due to 
interference between the Eq and multipoles, which 
vary rapidly at this energy due to the IV(1650)l/2 _ and 
_/V(1720)3/2 + resonances. Together with the slowly vary¬ 
ing E^ multipole, these three multipoles alone describe 
the E asymmetry quite well without the need for an addi¬ 
tional narrow resonance near 1.68 GeV. A similar analysis 
of the multipole content for the cos 9 cm bin centered at 
—0.6 shows that the interference of the Eq and M 2 mul¬ 
tipoles (the latter containing the V(1675)5/2 _ ) is respon¬ 
sible for the dip, with E+, E 2 and M 2 necessary to better 
approximate the full fit. Combined with the hints seen in 
Refs. [7J[8], the data presented here further motivate ad¬ 
ditional experimental investigations looking at other spin 
observables. 

In summary, we have presented the first measurements 
of the helicity asymmetry E in photoproduction of ry mesons 
from the proton. Initial investigation of these results with 
the Jiilich-Bonn dynamical coupled-channel approach show 
pronounced changes in the description of this variable when 
these new data are included, and demonstrate how these 
measurements are particularly impactful in constraining 
analyses of the excitation spectrum of the proton. With 
respect to the existence of an N^ + resonance near W ~ 1.7 
GeV suggested previously II El E|, the data obtained 
here do not demand the presence of such a state, but fur¬ 
ther measurements of other polarization observables would 
be helpful in gaining additional insight on that question. 
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Figure 2: Helicity asymmetry E for 7 p —¥ pp as a function of W 
at various values of cos$cm compared to several phenomenological 
predictions. 
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Figure 3: The helicity asymmetry E for 7 p —>■ pp using smaller W 
bins to explore the behavior of the helicity observable near W 1.7 
GeV. Predictions of the Jiilich-Bonn model as discussed in the text 
are shown by the solid line. 
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